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Abstract 

We study coefficients b n that are expressible as sums over the Li/Keiper 
constants Xj. We present a number of relations for and representations of b n . 
These include the expression of b n as a sum over nontrivial zeros of the Riemann 
zeta function, as well as integral representations. Conditional on the Riemann 
hypothesis, we provide the asymptotic form of b n ~ 2 -n ~ 2 Inn. 
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Introduction 

Let C denote the Riemann zeta function, and £(s) = (s/2)(s — l)7r~ s / 2 r(s/2)£(s) 
the classical completed zeta function, where T is the Gamma function [9l [T2l [TBI [T7] . 
Within the critical strip S, < Re s < 1, the complex zeros of ( and £ coincide, and 
we denote them by p. The ^-function is entire, of order 1, and of maximal type. 

Herein, we mainly investigate certain sums b n over complex zeta function zeros. 
We provide various representations and properties of these sums. We also supply 
some remarks on the Li criterion [TJ] for the Riemann hypothesis (RH). 

We recall that the Li equivalence for the RH results as a necessary and sufficient 
condition that the logarithmic derivative of the function £[1/(1 — z)] be analytic in 
the unit disk. This obtains from a conformal map of the critical strip to this disk. 
The equivalence [H] states that a necessary and sufficient condition for the nontrivial 
zeros of the ^-function to lie on the critical line Re s = 1/2 is that constants {A/ c }^ 1 
are nonnegative for every integer k. The sequence {X n }^ =l can be defined by 



(1.21 



The Aj's are connected to sums over the nontrivial zeros of ((s) by way of [TBI [T 

-H 

For further discussion of the Li criterion, its application, and results on series expan 
sion of the £ function, see for instance [21 HI EJ EJ [3j [7] . 
In particular, consider the expansions [TT] 



1 \ 00 X 00 

\ncp(z) = ln£ ( = - In 2 + £ -^z n = b + £ b n (z + 1)". (1.3) 



The middle expansion in terms of the Li/Keiper constants X n holds for \z\ < 5± < 1, 
where 5i = 1 corresponds to the Riemann hypothesis. Similarly, the right-most ex- 
pansion holds for \z + 1| < 82 < 2. Thus, (1.3) has overlapping domains of expansion, 
allowing analytic continuation. 

In [11], the coefficients b n of (1.3) are simply an inessential device. However, we 
treat their properties, including a "curious identity" 62 — ^3; in the next section. 
The latter relation is simply the beginning of an infinite set of relations that we 
make explicit. In the Appendix, we record approximate numerical values for the 
early coefficients. Reference [11] contains physical discussion, including attempting 
to regard £ quantum mechanical wave function. 

To emphasize that the 6 n 's are not required for the purposes of [TT], we have the 
following argument. Let N(T) be the count of nontrivial zeta zeros for < Im p < T. 
We have N(T) = 7r _1 Imm£(l/2 + iT), and it is well known that as T — > 00, N(T) = 

(T/2w) ln(T/27r) - T/(2vr) + O(lnT). Suppose that the RH holds. Then we have 

00 

A n = 2^(1 -cos ntij) > 0, (1.4) 
3=1 

where p = 1/2 ± ip,j or pj = (1/2) cot(#j/2). We can rewrite (1.4) as 

X n = 2 [1 - cos 6{p)]dN{p), (1.5) 
Jo 

where the lower limit just as well may be taken as p,\. Integrating by parts, we have 

roc AQ 

X n = -2n smn6—N(p)dp + 2(1- cosn6)N{p)\™, (1.6) 
jo dp 

where p{6) = |cot(6 l /2). For p — > 00, we have 6 = l/p + 0(l//i 3 ), 1 — cosn9 = 
—n 2 /2p 2 + 0(l/p 3 ), and the required limit on the right side of (1.6) is zero. We 
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therefore obtain the equivalent exact forms 

i-oo dQ 

X n = -2n sinn9—N(fi)dfi (1.7) 
Jo dfj, 

= 2n * smneNU)d6U). (1.8) 
Jo 



We can note that in (1.8), n/2 may be replaced by 2cot~ 1 (2/ii) ~ l/\/210, and 
therefore the range of integration is a relatively narrow one. 

We may quickly rewrite (1.7), as d8/dfi = — 4/(4/i 2 + 1). Furthermore, sinn# = 
sin#i7 n _i(cos#), where Uk is the kth Chebyshev polynomial of the second kind, and 
cos 6 = (A/j, 2 - l)/(4/i 2 + 1). We obtain 

thereby recovering (3.13) of [TT]. As observed there, on the RH, the values A n are 
indeed nonnegative. 

Implicit in [11] (p. 8) are the relations 

oo 

-ln2 = & + £&n, (1.10) 

n=l 

and for n > 1, 

^. . . AAA, 

l.ll 



J^n W J 



That (1.11) holds may be easily verified by using [TT] (3.4) and the orthogonality 
relation 

E(-irfrlf£l = (-!)'**», (i.i2) 
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with 5jk the Kronecker symbol. We shall have recourse to these relations in the 
following developments. 

Relations and representations of b n 

We have 

Proposition 1. For n > 1 we have 

1 



nz pes 



1-1 



l-2p 



(2-1) 



Corollary 1. We have &i = 0. 

Corollary 2. We have C'(l/ 2 )/C(V 2 ) = |(7 + f + 31n2 + lnvr). 

In (2.1) the sum includes zeros p along with 1 — p. (Owing to the functional 
equation of the £ function or ^-functions.) We write J2 P when the companion zero 
1 — p is explicitly taken into account. 

We have 



Corollary 3. 



Thus 



bn = -— E r ) E / 1 n- - (2-2) 



n2 n fr{ \jj ^ (1 - 2p)> 



[n/2] / x 



6n " n2»-i S l2JfcJ?(l-2p)2*- (2 ' 3) 



fc=2 V'V p 



We have 

Proposition 2. The coefficient b 2n+ i is always expressible as a rational linear com- 
bination of b 2n , b 2n -2, &2- 

Examples. We have 6 3 = 6 2 , ^5 = 26 4 — 6 3 , and 6 7 = 35 3 — 56 4 + 36 6 . 
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Let E 2fc = E pe .s (i-l p y^ • We have 
Corollary 4. We have the relation for n even 

1 ra/2-1 

— — F 

12" 



^n+l — 77^71 + 7 — 

2 n 



n 



n 

2k-l. 



(2.4) 



(2.5) 



k=l I 

We have 

Proposition 3. We have the summation relation for n > 2, 
Corollary 5. In particular, we have 62 = &3- 

Let be the Laguerre polynomial of degree n and parameter a pQ. Then we 
have the following representation. 
Proposition 4. We have 

Corollary 6. On the RH, with p = 1/2 + itj, and t,- is real, we have 

1 00 POO 

Write [H] 



(2.6) 



(2.7) 



^) = i+x;%^ (2.8) 

with £(1/2) = 1 + X^Li(— = ex P(M on the RH. The rapid asymptotic growth 
of cij with j has been described in |8J . We have 
Proposition 5. We have the recurrence relation for m > 1 

£(l/2)(m + l)& m = £(-1)^^ 



=ij= 



(2.9) 



Proposition 6. On the RH, we have for n > 1, 




Jo 



■*/2 sin[(n - 1)9/2] 
i cos™ +1 (#/2) 



N{y)d6. 



Corollary 7. On the RH, we have for n 3> 1, 

fe n ~ 2-"- 2 [In (n - 1) + 7 - 1 - 1h(4tt)] . 

Proposition 7. On the RH, we have 

\ n = — In n + (7 — 1 — In 2n)n + o(n). 
2 

A Corollary of Proposition 7 is Corollary 7. 

In the next section, proofs are supplied, as well as some discussion. 

Proofs of Propositions 

Proposition 1. We substitute the sum (1.2) into (1.11), 



The interchange of sums is justified by the absolute convergence of (1.2). 
Corollary 1 immediately follows as we have 






From the Hadamard product for the ^-function, we have 



m _ 1 



E 



1 



1-p/z 




Therefore, b\ = —\j (|), implying Corollary 2. 



Corollary 3 (2.2) follows by binomial expansion in (2.1). 

Remarks. Corollary 2 recovers what otherwise may be found by applying the 
functional equation of the ( function. 

Indeed, all odd order derivatives of £ are zero at 1/2. 

Proposition 2. We use (2.3) and put S 2fc = J2 P es (i-2p) 2fc • Then, for each n, S 2n 
may be eliminated between b 2n and b 2n +i, and the result follows. 
Corollary 4- We have for n even from Corollary 3 




(3.4a) 



and 





Therefore, from (3.4a) we have 




(3.5) 



We insert this equation into (3.4b) written in the form 




(3.6) 



We find 




Finally, we use a recursion r< 



relation for the binomial coefficient, ( m + 1 ) = (™) + ( n ™J, 



and obtain (2.4). 
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Examples. We have £ 4 = 8(-46 4 +36 3 ), h = 2b±-b 3 , E 6 = -12(256 3 -40& 4 + 166 6 ), 
and b 7 = (l/56)(-576 3 + 806 4 + 246 6 ). Of course, b 3 = b 2 = -(1/4)S 2 . 

Proposition 3. The result is a consequence of the functional equation of the £ 
function, so that 

=£(W) = (3-8) 

We have 



1-^7 U-l 



using 



Then 



oo oo / k \ 

&o + J>«(l + *) n E (- 1 )" _\ (l + ^) fc - n+1 , (3.9) 

n=l fc=n-l \ n V 



E (,\V + *)*- i+1 - (3-io) 



9 J fc f-ii v - x / 



hi <^)^° + S\l ( - ir C-i) (1+2)t 

= &0 + EE&n(-ir( n _ 1 ](l + ^. (3.11) 



fe=l n=l 

Comparing with the expansion (1.3), we obtain 



3=2 \3 ~ 1 

where we have used b\ = 0. Relation (2.5) follows. 
Proposition 4- We use (1.11) and [5] (Prop. 1) 



6 B = E(-l) i [" fW (3-12) 



A i = E / e-" u Lj_ 1 (u)du, (3.13) 



o 
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CO 



E / ^(-^ n E^(i) L U( u ) du - ( 3 - 14 ) 



so that 



We use the representation [lj (p. 286) with Bessel function J a for a > — 1 

n\ Jo 

giving 



L"(x) = / r+^JMVrt^dt, (3.15) 

n! io 



= e n u -i/2tl)l [°° e- 2t t n - l l 2 ,h(2^R)dt 

TV. JO 

= 2e u u- 1 ' 2 ^- / e~ 2a: x 2n J l (2^x)dx. (3.16) 
n! jo 

The integral is first evaluated [10] (p. 716) in terms of the confluent hypergeometric 
function 



7 \n. 



-l) n / U 



2«+i V 2 
(-l)M rl /« 



ii-i (|) • (3-17) 



2 n+l „ 

Here we have used Kummer's first transformation for the function iF% [1] (p. 191) as 
well as the relation 

L a n {z) = l^ + n °^ 1 F l {-n i a+\;z). (3.18) 
The insertion of (3.17) into (3.14) gives the Proposition. 
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Remarks. By integrating by parts we may verify that (2.6) returns the sum rep- 
resentation of Proposition 1. We have 

1 roc 

= "^gf e^A in(l , )A , (3 . 19 ) 

The integral converges since necessarily Re p > 0. We use the Laplace transform of 
a Laguerre polynomial [10J (p. 844), and we recover (2.1). 

The theory of Laguerre polynomials is pervasive in formulating the Li criterion 

By multiply differentiating (3.17), we have a family of summations, 

I^O^- 1 ^ = ^ L --> (I) • (3 - 20) 

Otherwise, we may follow the steps as above and find for a > — 1 and z ^ 1 



oc 



We have the contour integral representation 

„, s T(n + a + l) r(o+) ( z\ n . dt , 

where the contour encircles the origin in the positive direction and closes at Re z = 
— oo. This gives 

vt^'V M-tn 1 / (0+) (wgrl e tj (3 23) 

^ J ^ 2ttz 2^7-oc (l- M /2t)^ +1 t 2 ' 
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It may be verified that the residue at t — u/2 gives Ll l _ 1 (u/2)/n. 

The defining sums of Proposition 1 may be recovered from Corollary 6 in the 
following way Write on the RH 

d 



n2 

1 oo 

— Y 



1 f°° d 

— rV / e~ v cos(2tjv)— L n (v)dv 
n ~ l fyJo V 3 J dv K ' 



— e v cos(2tjv) 
dv 



L n (v)dv + lj . 

Then use [TU] (p. 846) for n even and odd to evaluate the integrals. 

Proposition 5. This follows from the identity — (p = <p'. We make use of 



(3.24) 



oo J 



i + e ctjZj = i + E «i E (i ) (-i)^^ + 1) 4 

j=i i=i «=o Vv 

oo oo oo /A 

= i + ec- 1 ^ + E E(-!r' ! + !)' 



(3.25) 



and similarly 



(p {z) 



E*V ; = EE(-1) H + 1)( J J% + i(^ + (3-26) 
j=i e=o j=e V7 

With some further series manipulations we obtain (2.9). 
Proposition 6. We have from (1.11) and (1.8) 

b n = 2E(~ir n / sinjeN(fi)d9. (3.27) 

We then apply the sum E^J-l)^)^' = (-l) n a; n (l + x)^ 1 at z = exp(±i0). 
Then simple manipulations yield (2.10). 

Corollary 7. In (2.10) we change variable to x = (n — 1)0/2. At the leading order 



in n we obtain 



>— n— 1 



1 r°° sin x 



7T JO X 



hi 



n — 1 

47TX 



- 1 



(3.28) 
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Performing the integral gives (2.11). 

Remarks. It is evident that (2.10) includes the cases 61 = 0, 62 = &3, and 65 = 



The asymptotic result (2.11) is consistent with the right-most expansion in (1.3) 
having a radius of convergence at most 2. 

Simply from the sum representation (2.1) one may suspect an asymptotic form 
b n ~ £c-~- One could also estimate b n from 



where the contour encircles z = — 1. 

For one of the integrals in (3.28), we first have x a sin x dx = cos(7ra/2)r(o; + 1) 
for —2 < Re a < 0. Then performing logarithmic differentiation and taking a — > — 1 
we obtain 



Numerically from (2.1) as a sum over the first 10 5 nontrivial zeta zeros pj)] we 
find 61000 - 9.21 x 1(T 302 while from (2.11) we have 6 100 o - 9.22 x 1(T 302 . See figure 
1 for a semilog plot of the first 1000 values of b n versus n (with h\ omitted). These 
numerical values, suggesting that indeed the right-most expansion in (1.3) has radius 
of convergence 2, could be taken as evidence that the RH holds. 

Proposition 7. Method 1. Similarly to Corollary 7, we have using (1.8), 



26 4 - 63- 




(3.29) 




(3.30) 




(3.31) 
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Therefore, at the leading order we have 



A„ 



n 



7T Jo 



smx 



x 



In 



n 
2ttx 



- 1 



dx. 



(3.32) 



with the error incurred being o(n). Using (3.30) gives the Proposition. 

Method 2. We alternatively use the expression (1.9) and have the expansions 



i 



i 



fi 12/x 3 

and sm(n9) = sin(n//x) + 0(l//i 3 ). We have 
sin(n0) (A/j, 2 + 1) . , 



O 



C/ n _i(cosi 
We then have 



sin 9 



4/i 



smmfc 



fi + O 



smmfc 



(3.33) 



(3.34) 



A„ ~ 32n 



30 H 2 N(ii) 
, (4/i 2 + l) 2 

= 2 

7T JO 



sin 



— ] d[i ~ 2n / 



sm — \ dfi 



Ml ^ 



N 

V 





sinf 



In 



sin vdv 
2nv\ 

+ 1 



n 



dv. 



(3.35) 



Using (3.30) we find (2.12). 

Corollary 7. We now reprove this Corollary as a result of Proposition 7. We have 
by (1.11), 



E(-ir 



,ri/ 



In j + 7 - 1 - ln(27r) + o(l) 



(3.36) 



In order to accurately approximate the summand, we use the digamma function 
ip = T'/r, and have = In j + o(l) for j ^> 1. Then we have, using an integral 
representation for ip [TD] (p. 943), 

bn = E(- 1 ) J ~ n Q + 7] + I - 1 - M27T) + 0(1)} 
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= £(-ir n 

j=n 

1 r 1 



3 



2 Jo 



t 



n-l 



1 r l (P- 1 - 1 

2 Jo \ t- 1 



dt+ | - 1 - 1h(2tt) + o{l) 



(t + l)n+l 2 ™ +1 



1 1 



2 2 n+1 



^(n) + 7 — In 2 — 



(f-1) 
1 



+ 



7 



1 - ln(2vr) 



n 



+ p - 1 - ln(27r)] 

l 2 On+l 



2n+l 2 n+1 

1 o(l) 



^lnn + + ( 7 -l-ln7r 



In 2 



2«+i 1 2 n+1 
1 ofl) 



2n+l 2 n+1 



(3.37) 



Remarks. The result (2.12) is not new [5], but we include it and the method of 
proof as a companion to Proposition 6 and Corollary 7. We suspect that the o{n) 
terms in (2.12) are of size 0{n l ' 2+e ) for any e > 0. 

Generally alternating binomial sums may be difficult to estimate, but we have 
done so in recovering Corollary 7. 

Regarding (3.36) and (3.37), it is possible to use an even more accurate approxi- 
mation to In j, with ^(J + l/2) = In j + 0(1/ j 2 ), but at the cost of a more complicated 
integral to perform. 

Suppose that Aj has a subdominant term close to yfj. Then we expect there to 
be a correction term in b n close to 



EC- 1 ) 



-n 3 



j-n 



\ n J 3 



1/2 



EC- 1 : 



j—n I 3 



]=n 



Tin) 



T(n+ 1/2) 



2 Fi [n, n + 1; n + -; — 1 ~ 



n) T(j + 1/2) 
\[2 1 



2 n+l 



?? 



(3.38) 



Here, iF\ is the Gauss hypergeometric function [H [10] , and by transformation rules 
[TO] (p. 1043), the 2-F\ function in (3.35) is the same as 



-n-l 



. 1 1 1 

a F 1 (n,-;n+-;- 



2- n ~ 1 V2 2 Fi 
15 



11 1 , > 

--, -: n H — ; — 1 . 

2 2 2 



(3.39) 



The above argument extends so that if Xj has a subdominant term 
expect in b n a term close to 
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Appendix: Values of b n 

Exact expressions for b n can be written in terms of ln7r, polygammic constants 
■0C ? ')(l/4), and the derivatives £^(1/2). The following table gives approximate nu- 
merical values for the initial 6„'s. 



n 


K 





-0.698922 


1 





2 


0.00144406 


3 


0.00144406 


4 


0.00108297 


5 


0.000721886 


6 


0.000451088 


7 


0.00027058 


8 


0.000157786 


9 


0.0000901269 


10 


0.0000506726 


11 


0.0000281364 


12 


0.0000154657 


13 


8.43018 x 10~ 6 


14 


4.56299 x 10~ 6 


15 


2.45502 x 10~ 6 


16 


1.31 x 10~ 6 


17 


6.99 x 10~ 7 


18 


3.71 x 10~ y 


19 


1.96 x 10~ 7 


20 


1.03 x 10~ 7 


21 


5.44 x 10~ 8 


22 


2.85 x 10~ 8 


23 


1.49 x 10~ 8 


24 


7.79 x 10- 9 


25 


4.06 x 10~ 9 



The values b , . . . , b 15 have been obtained in Mathematica by series expansion of 



17 



the if function, (1.3). The remaining values have been found in Matlab by summing 
over the first 10 5 complex zeta zeros [T5j, (2.1). 
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